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Abstract:  Root-based  direction  finding  algorithms  (DF)  have  several  advantages  over 
search-based  DF  algorithms.  A  key  advantage  is  the  fact  that  they  do  not  require  the  array 
steering  vector;  this  is  because  these  algorithms  presume  equalized  element  radiation  patterns. 
In  this  paper  the  WIPL-D  code  is  used  in  designing  an  array  of  rectangular  probe-fed  patch 
antennas  with  equalized  radiation  patterns  for  measuring  the  range  and  bearing  of  RF  emitters 
in  the  PCS  band  (1900-1920  MHz).  Direction  of  arrival  (DoA)  estimation  results  based  on 
simulations  and  measured  data  are  presented  and  used  as  a  measure  of  element  patterns 
deviation  from  equality. 


1.  INTRODUCTION 

Direction  finding  (DF)  algorithms  for  DoA  estimation  is  a  topic  that  has  been  studied 
thoroughly  in  the  past  few  decades  mainly  by  researchers  in  the  signal  processing  and  antenna 
theory  communities.  From  signal  processing  point  of  view  the  focus  has  been  on  maximizing 
the  number  of  DoAs  that  can  be  accurately  estimated  and  at  the  same  time  reducing  the 
computational  cost  involved  in  this  work.  Several  algorithms  exist  and  can  be  used  to  estimate 
the  DoA  of  incident  signal  on  an  antenna  structure  [l]-[3].  DF  algorithms  are  classified  as 
either  search-based  or  root-based.  We  refer  to  the  former  class  as  S-DF  algorithms  and  to  the 
later  class  as  R-DF  algorithms.  R-DF  algorithms  have  several  advantages  over  S-DF  algorithms. 
In  addition  to  their  less  computation  cost  where  DoAs  are  calculated  by  finding  the  roots  of  a 
polynomial  of  certain  order  rather  than  going  through  intensive  search  in  the  whole  angular 
domain,  a  key  advantage  from  using  R-DF  algorithms  is  the  fact  that  they  do  not  require  the 
array  steering  vector.  This  is  because  these  algorithms  presume  equalized  element  patterns  and 
their  accuracy  depends  on  how  much  the  element  patterns  deviate  from  equality. 

This  paper  investigates  the  design  and  performance  analysis  of  antenna  structures  with 
equalized  element  patterns  for  R-DF  algorithms.  Adding  a  number  of  passive  elements  around 
the  center  active  elements  and  terminating  them  using  a  suitable  set  of  loads  minimize  the 
deviation  of  the  element  patterns  of  the  center  elements  from  equality  [4].  Return  loss  and 
element  patterns  of  two  antenna  structures  comprised  of  a  number  of  rectangular  probe-fed 
patch  antennas  are  calculated  using  the  WIPF-D  code  [5].  The  first  structure  consists  of  four 
elements  and  no  passive  elements.  In  the  second  structure  three  passive  elements  are  added  on 
both  sides  of  the  four  active  elements  and  were  terminated  in  50Q.  DoA  estimation  results 
using  the  two  antenna  structures  are  compared  with  DoA  estimated  using  measured  data.  The 
DoA  accuracy  is  used  as  a  measure  of  how  much  element  patterns  deviate  from  equality.  The 
paper  is  organized  into  four  sections.  Following  this  introductory  section,  in  section  two  return 
loss  and  element  patterns  calculated  using  WIPF-D  are  presented  for  the  two  antenna  structures 
mentioned  before.  DoA  estimation  results  using  either  WIPE  calculated  patterns  or  measured 
data  are  compared  and  presented  in  section  three.  Finally  conclusions  are  provided  in  section 
four. 


2.  EQUALIZATION  OF  ELEMENT  PATTERNS 

As  shown  in  figure  (1),  equalization  of  element  patterns  is  done  through  adding  a 
number  of  passive  elements  around  the  middle  active  elements.  The  top  structure  in  figure  (1) 
is  for  an  array  of  4-element  of  patch  antennas  with  no  passive  elements.  In  the  bottom  structure 
3  passive  elements  on  each  side  of  the  active  elements  were  used,  respectively.  We  refer  to  the 
top  and  bottom  structures  as  Arrayl  and  Array2,  respectively.  Figure  (2)  (left)  shows  top  and 


©2005  ACES 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

01  JAN  2005 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Equalization  of  Numerically  Calculated  Element  Patterns  for  Root-Based 

5b.  GRANT  NUMBER 

uu  cuiuii  r  uiuiug  /iiguiiuiins 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Electrical  and  Computer  Engineering,  The  George 
Washington  University  Washington  DC  20052,  USA 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM001846,  Applied  Computational  Electromagnetics  Society  2005  Journal,  Newsletter,  and 
Conference.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

4 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


bottom  view  of  a  4-element  array  of  patch  antennas  modeled  in  WIPL.  The  dimensions  of  the 
patch,  ground  plane  and  coax  were  identical  to  those  given  in  [6].  The  inter-element  spacing  for 
Array  1  and  Array2  was  kept  fixed  and  equal  to  X/2.  The  bottom  and  top  views  of  the 
10-element  array  of  patches  (Array2)  are  shown  in  figure  (2)  (right).  The  whole  structure  is 
modeled  in  the  same  way  as  with  Arrayl  except  that  passive  elements  were  terminated  in  50Q 
loads. 

4-element  array 
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Fig.  1, 4  -element  array,  top,  and  10-element  array  with  3  passive  elements  on  each  side,  bottom 


Fig.  2  A  4-element  (left)  and  10-element  (right)  array  of  patch  antennas  modeled  in  WIPL 


The  return  loss  in  the  (1.9-1.92)  GHz  frequency  band  of  Arrayl  and  Array2  calculated  using  the 
WIPL-D  code  and  are  shown  in  figure  (3).  For  both  arrays  it  is  observed  that  the  resonant 
frequency  of  each  of  the  middle  elements  is  tuned  and  shifted  by  1MHz  compared  to  the  outer 
elements.  Element  patterns  of  Arrayl  calculated  at  1.91  GHz  are  presented  in  figure  (4),  for 
fixed  cp  and  fixed  0,  respectively.  Results  for  Array2  are  presented  in  figure  (5).  It  is  clear  from 
the  figures  that  the  element  patterns  of  the  four  elements  significantly  deviate  from  equality. 


Fig.  3  Return  loss  of  Arrayl  (left),  Array2  (right) 


Fig.  4  element  patterns  of  array  Arrayl  at  cp=7i/4  (left)  and  at  0=44°  (right) 


Fig.  5  Element  patterns  of  the  active  4  elements  of  Array2  at  cp=7i/4  and  at  0=44°  (right) 

The  accuracy  of  DoA  estimation  using  the  element  patterns  of  Array  1  and  Array  are  presented 
in  the  following  section  and  are  compared  to  DoA  results  based  on  measured  data  with  Array2. 


3.  NARROWBAND  ROOT-BASED  DIRECTION  FINDING  ALGORITHMS 

In  this  section  accuracy  of  DoA  estimation  using  the  modified  root-Pisarenko  (MRP) 
[7]  is  analyzed.  The  analysis  will  be  done  for  Array  1  and  Array2  and  are  compared  to  measured 
data  collected  from  an  antenna  system  consisting  of  Array2. 


Fig.  6,  Geometry  of  the  antenna  array 


Referring  to  figure  (6),  the  narrowband  signal  model  is  assumed  and  is  described  with 

z(t)  =  As(0  +  v(^)  (1) 


Where  z(t)  is  the  output  of  the  I-Q  channels,  A  is  the  steering  vector,  s  is  the  unknown  source 
signal  and  v  is  vector  of  additive  Gaussian  noise.  The  number  of  antenna  elements  is  N  and  the 
number  of  sources  is  L.  Elements  of  the  steering  matrix  are 


wt-“-  (2) 


In  equation  (2)  fn  is  n-th  element  pattern  in  the  terminated  array  environment  including  mutual 
coupling  effects.  For  uniformly  spaced  linear  arrays  along  the  x-axis,  the  steering  vector  is 

Where  ^  =  sin  9  cos  cp  and  d  is  the  inter-element  spacing.  When  element  patterns  are  equalized 
the  steering  vector  in  equation  (3)  reduces  to 

«  M  .  -  J_  /  M  [l  ^  ...  JW-'Wf  (4) 

Table  1  summarizes  results  of  DoA  estimation  using  MRP  algorithm  and  Uniform  Linear  Array 
(ULA)  consisted  of  4  active  and  6  passive  rectangular  probe-fed  patch  antennas.  The  true 
position  measured  by  transit  is  19.18°.  In  DoA  estimation  from  measured  data  one  ULA  with  6 
passive  and  4  active  elements  was  used.  The  array  was  a  part  of  the  RF  emitter  range  estimation 


system  comprised  of  two  antenna  arrays.  In  Do  A  estimation  from  simulated  data  Array  1  and 
Array2  configuration  were  used  with  element  radiation  patterns  calculated  numerically  using 
the  WIPL-D  code  [5].  ULA  configuration  with  passive  elements  exhibited  better  Do  A  accuracy 
in  both  cases  with  measured  and  with  simulated  data. 


Fig.  7  DoA  estimation  roots  using  Array  1  (left)  andArray2  (right) 


Estimated  from  measured  data  using  Array2 

19.86° 

Estimated  from  simulation  using  Array2 

20.25° 

Estimated  from  simulation  using  Array  1 

21.61° 

Table  1  DoA  estimated  from  simulations  and  measured  data  using  Array  1  and  Array2 


4.  CONCLUDING  REMARKS 

An  array  of  rectangular  probe-fed  patch  antennas  with  equal  radiation  patterns  for  root-based 
direction  finding  algorithms  was  modeled  using  the  WIPL-D  code.  Adding  three  passive 
elements  on  both  sides  of  four  active  elements  and  terminating  them  in  50Q  sufficiently 
equalize  patterns  of  the  active  elements.  DoA  estimation  results  using  this  array  with  equalized 
patterns  shows  the  accuracy  improvement  when  compared  to  array  without  passive  elements. 
Finally,  the  DoA  estimation  results  obtained  from  simulated  data  agree  with  DoA  estimation 
results  obtained  from  measured  data. 
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